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Abstract 

In the present work, an exact solution of Einstein's field equations 
representing the field of a binary system, is obtained. It is shown 
that this solution is a time-dependent one. This solution is a modified 
version of that obtained by Curzon, since it reduces to the later one in 
the static case. The equations of motion for a test particle in the field 
of a binary system are formulated and solved. Such equations can be 
used to study the motion of a third body, or a photon, in the field of 
a binary system. 

Introduction 

What is meant by the two body problem is the problem of two structurless 
non-spinning point-like particles, characterized by two mass parameters mi 
and m 2 , moving under their mutual gravitational interaction. In order to 
formulate this problem, one aims to get an explicit expression of accelera- 
tion of the binaries in terms of their positions and velocities. To discuss this 
problem we can dissolve it into two aspects: i) obtaining the equations of 
motion of the two interacting bodies, ii) solving these equations. 
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In the context of Newtonian gravity the members of the system are consid- 
ered as widely separated objects, such that the contribution of the non-linear 
effects can be neglected. Newtonian gravity has a linear structure that en- 
ables us to derive the equations of motion of the binary components, and also 
to get an exact solution for these equations. It gives a full treatment of the 
binary system. On the other hand, in General Relativity, The field equations 
have a non-linear hyperbolic structure, so that it is not easy even to derive 
the equations of motion for such systems. Since in General Relativity the 
equations of motion are embedded in the field equations, consequently it is 
very difficult to derive the field equations as a linear functional of the matter 
distribution independently of the equations of motion [TJ [5]. 

One may ask why we are obliged to deal with General Relativity. Actually 
the Newtonian theory is very limited, since it deals only with systems, have 
a widely separated and slowly moving objects. Therefore, this theory can 
not predict the behavior of systems with a small separation and fast moving 
objects. Moreover it can not account for the accumulating small feed back 
effects, e.g. the advance of perihelia of planetary orbits. Fortunately, in such 
cases we are saved by General Relativity, which gives a good agreement with 
the observational and experimental tests. In order to avoid the above men- 
tioned difficulties, most of the relativistic treatments in the literature have 
used an approximation technique, i.e. the Parameterized Post Newtonian 
(PPN) method. This method expands the relativistic effects in power series 
of v/c, when only the relativistic expressions with order of v 2 /c 2 added to 
the Newton's law, so it is called first post-Newtonian (1 PN) approximation 
|lUj . Recently, the equations of motion of binary systems have been derived 
at the third and half post-Newtonian (3.5 PN) order. However the beauty of 
this technique is its applicability for various classes of relativistic theories of 
gravity and also on its attractive mathematical formalism. The weak point 
of this technique is the breaking of the general covariance principle. 

Exact solutions in General Relativity have played an important role to 
give a full description of some physical problems (e.g. Schwarzschild, kerr, 
Reissner-Nordstrom) . The meaning of an "Exact Solution" is that the metric 
is represented by a coordinate system, in terms of the well-known analytical 
functions. Unfortunately, there are few exact solutions for real physical prob- 
lems that have been established, as mentioned by Kinnersky c.f. [B] "Most of 
the known exact solutions describe situations which are frankly unphysical" . 
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Although we have a lot of these non-physical solutions, one of the most 
important problems, which remained without an exact solution, is the phys- 
ical field of two bodies. In general, for any physical problem we aim to 
construct a mathematical model; using some reasonable conditions; defined 
by a certain set of differential equations. Actually it is not easy to interpret 
the obtained solution in some physical theories such as General Relativity, 
because of its high non-linearity. But even if we could not understand the 
qualitative features of an exact solution, we can compare it with approximate 
results, which will be very useful, to check the validity of this solution. 

The aim of the present work is to find an exact solution, of Einstein's field 
equations, which can be used to represent the gravitational field of a binary 
systems. This certifies a difficult task, but we are going to simplify this prob- 
lem in the following manner. In the literature there are many trials to solve 
two body problem in the frame work of General Relativity. The solution 
required is expected to represent the field of two revolving objects (singular- 
ities) in one system. In 1924 Curzon obtained a solution, to Einstein's field 
equations, representing the field of two static singularities. Unfortunately, 
his solution has been classified among un-physical solutions, since it does not 
represent a real astrophysical system, we are going to adjust this solution 
in order to be suitable for describing a binary system. For this reason, we 
give a brief review of Curzon solution in §[TJ with its necessary background. 
In section we are going to modify (adjust) this solution to represent the 
field of a real binary system. We test the solution by reducing to the flat 
space and schwarzschild solutions, under certain conditions, in §|3} Motion 
in the modified Curzon field is treated in section §|4| The work is discussed 
and concluded in £j5j 

We assumed that the two members of this system have the same masses 
and rotate in a circular orbit about the center of mass. We assumed further, 
that the line of sight is in the plane of the orbit. This represents a two 
body problem which represents many physical configurations especially PSR 
J0737-3039 jSEj. 
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1 Curzon Solution 



This solution is a Weyl class one, which was found soon after the birth of 
GR. It refers to two "monopoles" on the axis of symmetry. This is mentioned 
by Bonnor |2J as "Probably the most perspicacious of all exact solutions in 
GR". 

1.1 Standard Curzon Static Field 

In this section, we aim to review the solution, obtained by Curzon (1924), of 
Einstein's equations 

R a p = 0, (1.1) 

for an axial symmetric gravitational field produced by two singularities on 
the axis x 1 separated by a distance 2a. If the origin of a reference frame 
lies at the mid point between the two singularities on this axis, so for an 
arbitrary point P in the plane of x 1 and x 2 coordinates the direction of these 
singularities can be defined using bipolar coordinates f\ and f2, as shown in 
Figure 1: 




Figure 1: The arbitrary point P on the plane of X\ and X2 co- 
ordinates is referred to the singularities m x and m 2 by bipolar 
coordinates fx and r 2 , respectively. 
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where 



f\= (x 1 -a) 2 +{x 2 ) 



rl=(x l + ay+(x 2 )' 
The cylindrical coordinates used in this study are 



x 1 



x 2 



P, 



x 3 



x 4 



t 



(1.2) 



(1.3) 



The metric characterizing the space with axial-symmetric static gravitational 
field H]. 

ds 2 = -e* (dz 2 + dp 2 ) - e~ p p 2 d<f) 2 + e^dt 2 (1.4) 

where /2 = p,(z, p) , 
9 = u(z,p). 
So that the metric coefficients g a p are 



9n = ~e h 



-2 -P 

#33 = -p e , 



#44 



922 = -er 

9a(B = 0, 

The components of the contravariant tensor g af5 are given by 
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Using the definition of Christoffel symbol 

def 
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we get for the space (|1.4l) the following non- vanishing components, 
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(1.5) 



i - , du _ . dp _ dv dv 

where p\ denotes — , p 2 denotes — , v\ denotes — and v 2 denotes — . 
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dx a dx 13 
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dp' 
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Now by using Einstein's field equation in free space (1.1), 
«/■ d 2 log 



d 



dx'r 



dx~i [a/3 



7 



7 
ae 



e 

7/3 



= 0, 
(1.6) 



the following set of differential equations is obtained 



Rn 


=\{fln + P>22 + v\ + P2IP) = 0, 


(1.7.1) 


R\2 


= §0^2 - V X /p - pi/p) = 0, 


(1.7.2) 


R22 


=|(a*ii + A*22 + v\ - p 2 /p - 2v 2 /p) = 0, 


(1.7.3) 


R33 


= - |e-^-y (z/ n + v 22 + z/ 2 /p) = 0, 


(1.7.4) 


-R44 


= - |e^ +p (i7 n + z> 22 + iAj/p) = 0. 


(1.7.5) 



Equation (1.7.5) in this set is a mere repetition of equation (1.7.4). Thus 
we have only four equations for p and v. Since v is a function of z and p, 



equation (1.7.5) can be written in the form, 

V 2 v = 0, 

where V = £p + 1 f = + ^2 i n the cylindrical coordinates symmetrical about 
the z-axis. It is clear that the above equation represents the classical (flat- 
space) Laplace's equation, hence the solution of this equation is 



v 



-2 



mi m 2 
f x f 2 



;i.8) 



where mx and m 2 are constants of the integration. One may realize how the 
contribution of the classical theory of gravity is represented by this linear 
equation. To obtain the relativistic effects of the curved space-time, let us 
introduce a new function A, such that 



A = p + v, 



1.9) 



by using (1.7.2), we get 



Al = px + Vx = V\V2p, 



1.10) 
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by using (1.7.1) and (1.7.3), we get 

A2 = A*2 + ^2 

Now, 



d\ = Xidz + X 2 dp, 
hence, substituting from ( |1.10[ ) and (1.11) into (1.12), we get 



(1.11) 
(1.12) 



,2 ( - 



dz + 2p [p 2 



+ (4p 2 (z + a)dz + 2p [p 2 



- a) 2 } dp) ^ 
' 1 

[z + a) 2 } dp) ^ 



+4 (z p dz - [ 



z 2 



p 2 ] dp) mim2 



3^3 



;i.i3) 



1' 2 



Since the metric of the manifold given by equation ( |1.4 ), thus the integration 
of equation (1.13) should reduce to the flat space at infinity as f\ and f 2 — > 



00, so A tends to zero. 



A 



1 1 



m 



2 7,2 



mim 2 
a 2 



z 2 + p 2 



r x r 2 



1 



;i.l4) 



One can realize that the function A is determined from v by quadrature up 
to an additive constant. Also, the metric function A has the property that 
lim A 7^ 0, so there is some kind of a compression described by a canonical 

singularity between the particles holding particles apart [9] . From equations 
(1.8), (1.9) and (1.14), the function p is given by 

777.17712 



777i 7772 
fi f 2 



1 2 
777| 7772 

-4 -4 



P 2 + 



z 2 + p 2 



r x r 2 



;i.i5) 



Now we define two angles ol\ and a 2 , as the angles between by AP, BP and 
the axis z, respectively (see Figure 1). Finally the solution of (1.1), in the 
present case, can be written as 



777! 777 2 

f\ f 2 



1.16) 



mi m 2 
fx f 2 



2 2 
TTli 2 _ ^^2 2 — 

sin ax + -jz- sin a 2 



r 



n 



mxm 2 . 2 (a 1 -a 2 \ 
-2 — S777 { 2 ) ■ 



;i.i7) 
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The solution given by (1.16) and (1.17) is the solution given by Curzon jl] 
in its standard form. 



1.2 Stability of Curzon Solution 



The relativistic solution, given above, can be rewritten in the following equiv- 
alent form [9]. 



e^dt 2 - e 



where 



mi m 2 
fx f 2 



e~ x (dz 2 + dp 2 ) + p 2 



A 



m\ _ 2 m 2 _ 2 171x1712 
^4~P -4 P 



r. 



z 2 + p 2 — a 2 
fxf 2 



1 



To compute the force between the two singularities, one easily see that A(0) 
along the z-axis between the two singularities, can be given as 



A 



mxrri2 



as 



171x1712 



then the stress force between the two masses 

GMxM 2 
(2D) 2 ' 



;i.l8) 



1.19) 



if we take the Newtonian limit D ^> mx,rri2, the medium between the two 
masses contains a compression merely the Newtonian force attraction [9]. 



2 Modification of Curzon Solution 

Actually the static case, which had been studied by Curzon, does not rep- 
resent any physical configuration (e.g. real binary systems in Nature) in 
view of the fact that the binary systems are always dynamical systems. The 
main aim of the present section is to show how to modify Curzon solution 
such that it can be used to represent the field of a real physical binary system. 
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In order to describe such systems, the effect of time should be introduced, 
by considering a co-rotating coordinate system, associated with the angular 
motion of the collinear singularities, with respect to a frame of reference fixed 
at the center. It is similar to the situation of using co-moving coordinate 
systems in cosmological applications. 

2.1 General Outlines of the Modification 

Curzon has chosen the point P in the (z-p) plane. Since A and B are static 
they could not define a particular plane, he always can reorient the coor- 
dinates by choosing a particular value of an angle to make ABP plane 
always coincides with the (z, p) plane, without affecting generality. While 
the non-static case, such collinear singularities A and B are rotating around 
a common center C by an angular velocity u. It is clear that the moving 
singularities are defining a particular plane (orbital plane), so the point P 
should be an arbitrary point, not necessary in this plane. The separations 
between this point and the two singularities (A, B) are respectively 



where the angle represents the inclination of the orbital plane to the plane 
of the sky. As (f) — > 7r/2 the binary system tends to be an eclipsing binary 
system with respect to an observer at the point P. 

In general the singularities (binary components) are assumed to be sepa- 
rated from C by distances a and b (a ^ b), and the angular velocity u in this 
case is a function of time. To simplify the problem we are going to assume 
that the motion is circular as following: 

1. The masses of two singularities are comparable masses. 

2. The distances, a and 6, are equal. 

3. The angular velocity oj of the collinear singularities is constant. 

4. The orbits of the two singularities will coincident on each other pro- 
ducing a circular orbit with a radius a. 




(2.1) 
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Now before discussing the non-static case, we will return to standard form 
of Curzon metric (1.4), and show its form in other coordinate systems. This 
is done in order to facilitate comparison with some special cases. 



2.2 A Particular Choice of Coordinate Systems 

In order to refer the points of the manifold to the Cartesian coordinate sys- 
tem, as general covariance is still preserved since we use tensors in the for- 
malism, we use the transformation 

TI: (z,p,(j),t) — > (x,y,z,t) 



i.e. 



z — z, 

p = a/x 2 + y 2 , 
4> = tan _1 (x/^), 
i=t. 



> 



(2.2) 



Applying (2.2) to Curzon metric (1.4) and recalling that ds is a scalar, we 
get 



ds 2 (x p ) = ds 2 (x a " 



So we can rewrite (1.4) in the form 

x 2 e^ + y 2 2 



ds 2 



x 2 + y 2 

y 2 e^ + x 2 e~ v 
x 2 + y 2 



r/.H - — — ^-^{e' 1 — e u ) dx dy 
x l + y 2 

dy 2 - e A dz 2 + e p di 2 



where 



(2.3) 



v[x,y,z) 



-2 



mi m 2 
fi f 2 



(2.4.1) 



p{x,y,z) = 2 


mi 


m 2 
f 2 _ 


- {x 2 + y) 


mi m 2 

~4 ~A 




^x 


2 + y 2 + z 2 - 




' a 2 




fif 2 


°1 -1 



(2.4.2) 
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and 



f\= x 2 + y 2 + z 2 + a 2 - 2az 
r\= x 2 + y 2 + z 2 + a 2 + 2az 



(2.5) 



The metric (2.3) represents the gravitational field of two singularities, and 



has axial symmetry about the z-axis. 

Again we are going to apply a second transformation to represent the 



metric (2.3) in spherical polar coordinate, as follows 



TIT. (x,y,z,i) 



(r,9,(f>,t) 



i.e. 



> 



(2.6) 



x = f cos 9 
y = r sin 9 cos < 
z = f sin 9 sin ^ 
t = i 

where < 9 < tt and < <fi < 2n. One can notice that the transformation 
Til is not the conventional transformation. Consequently, the angle <fi is not 
the same angle 0, given in the standard Curzon solution. 
Since, 

ds 2 (x~<) = ds 2 (Z a ). 



Then, the metric coefficient will be given by 



9n 



922 



923 



933 



914 



~ cos 2 9 sin 2 (j> + cos 2 6 e 

r 7, 7. 7. 

cos 2 9 + sin 2 9 cos 2 <f> 
f 2 sin 9 sin 6 cos 9 cos 



cos 2 9 + sin 2 9 cos 2 - 



—f sin 9 



cos 2 9 sin 



+ cos 



cos 2 9 + sin 2 9 cos 2 ■ 



(2.7.1) 
(2.7.2) 

(2.7.3) 

(2.7.4) 
(2.7.5) 



11 



where 



u(f,9,- 



p.(r,6,<j>) 



-2 



mi ^2 



(2.8.1) 



mi m2 



— f 2 ^cos 2 9 + sin 2 cos 2 

171x1712 



mi m 2 



a' 



f 2 — a 2 



rxr 2 



- 1 



and 



(f + a) 
(r + a) 



2af(l + sin# sin< 
2af(l — sin# sin< 



(2.8.2) 



(2.9) 



In view of the previous illustration, considering the non-static case assump- 



tions, given in ^21, we apply the third transformation 

Till: (r,§,lt)^(r, 0,0,*) 



i.e. 



r = r, 



(2.10) 



> 

<p = 4> + u;i, 
* = *, 

which if combined with the transformation law of the scalar, 

ds 2 {x a ) = ds 2 {x^), 

would give the following non-vanishing components of the metric tensor in 
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terms of the new coordinate system (r, 9, 0, t). 

2 cos 2 9 sin 2 (0 + cut) + cos 2 (0 + cut) e 



922 



923 



924 



933 



934 



944 



[cos 2 9 + sin 2 9 cos 2 (0 + cut) 

r 2 sin 9 sin (0 + cut) cos 9 cos (0 + out) 
[cos 2 # + sin 2 9 cos 2 (0 + cut 



(2.11.1) 
(2.11.2) 

(e^-e-"), (2.11.3) 



— cu 



r 2 sin # sin (0 + cut) cos 9 cos (0 + cut) 
[cos 2 9 + sin 2 # cos 2 (0 + cut 



9 9 „ cos 2 6 1 sin 2 (0 + out) e~ u + cos 2 (0 + cut) 
-r sin^ — 



(e^-e^) ,(2.11.4) 
(2.11.5) 



[cos 2 9 + sin 2 9 cos 2 (0 + cut)] 

2 . 2 „ cos 2 6> sin 2 (0 + cut) e _i/ + cos 2 (0 + cut) e M . , , . 
our sm 9 F — — —ftz — - , (2.11.6) 



[cos 2 9 + sin 2 9 cos 2 (0 + cut 

- [cuV sin 2 9 (cos 2 9 sin 2 (0 + cut) e~ v + cos 2 (0 + cut) e M ) 

- (cos 2 9 + sin 2 9 cos 2 (0 + cut)) e v ] / 

[cos 2 9 + sin 2 cos 2 (0 + cut)], (2.11.7) 



while (|2.8.1|), (|2.8.2|) and ([2.9} now read: 

-2 



u(r,6,(/),t) 



fi(r,9,(p,t) 



mi m 2 
ri r 2 



(2.12.1) 



mi m 2 
n r 2 



-r 2 (cos 2 9 + sin 2 # cos 2 (0 + cut)) 
mim 2 



m 



_L |_ ^ 
4 4 



+ 



r 2 — a 2 



rir 2 



(2.12.2) 
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and 



(r + a) 2 — 2ar{\ + sin 6 sin (0 + wi)), 
(r + a) 2 — 2ar(l — sin 6 sin (0 + cot)). 



(2-13) 



Thus the set of quantities (2.11.1 )-(2.11.7) represents the field which is pro- 



duced by a dynamical two-body system, at any chosen point outside the two 
singularities. 

3 Check of the modified solution 

The line element of the modified Curzon solution is given in terms of spherical 
polar coordinates (r,8,(f),t). In order to check this solution one should obtain 
the flat space metric and Schwarzschild metric as limiting cases under certain 
conditions. 



3.1 First check: flat space metric 



It can be easily shown that the \i and v functions, given by (2.12.1) and 



(2.12.2), vanish as r 1; r 2 — > oo, i.e. at a large distance from the binary 



system. In this case the metric coefficients (2.11.1 )-(2.11.7) will become 

5-n = -1, 5-22 = ~r 2 , #33 = -^sin 2 ^, g 4i = 1 
and the line element will reduce to 

ds 2 = -dr 2 - r 2 d6 2 - r 2 sin 2 Odcf 2 + dt 2 
which describes the space in absence of the gravitational field. 

3.2 Second check: Schwarzschild metric 

Taking u = 0, mi + m 2 = m, and taking r large enoug h, r\ m r 2 ~ r, such 



that O 



— > 0. In this case the metric coefficients (2.11.1 )-(2.11.7) will 



become 

gn = -e M , g 2 2 = -rV, g 33 = -r 2 sin 2 9e^, g M = e u , 
and the line element can be written as 

ds i = - e 2 ^ (dr 2 + r 2 d6 2 + r 2 sin 2 6d<p 2 ) + e 2 »*dt 2 . 



14 



Where 



m 1 


r 

+ - 




n 


7 


i + 


m 2 






T2 



m 
r 



+ o 



then the line element can be rewritten as 



ds 2 







m 


+ o 






1 + 










r 










m 


+ o 




+ 


1 + 










r 







m 
r 



(dr 2 + r 2 d9 2 + r 2 sin 2 



dt 2 



Apply the transformation R = r + m, we get 



ds 2 



-i 



dR 2 -R 2 d6 2 -R 2 sin 2 



2 +(i-^^ 2 , 



i2 / 



which represents the line element of Schwarzschild gravitational field in its 
standard form. 



4 Motion in The Modified Curzon Field 

"spacetime tells matter how to move, and matter tells spacetime how to 
curve" [8j. We have seen how matter tells spacetime how to curve, now we 
would like to search how spacetime tells matter how to move!! So we calculate 
the non-vanishing Christoffel symbol coefficients of the second kind (sym- 



metric in first two indices) for the space represented by (2.11.1 )-(2.11.7). By 
using the calculated values of the Christoffel symbols and apply the geodesic 
equation, 

d 2 x a ^ dx^ dx v 

1- y = 0, 

ds 2 ^ v ds ds 

one can formulate the equations of motion of a test particle in the gravi- 
tational field of a binary system. In this way we can describe the motion 
of a test particle like a third body in the binary system, without adopting 
perturbation technique, or we can describe the motion of a massless particle 
(e.g. photon) in the field of the binary pulsar. To test planer motion, we 
put x 2 = 9 in the equations of motion and using the above values of the 
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Christoffel symbols, we can see that the differential equation for the angle 9 
can be written as, 



<P6 
ds- 



+ 



,dr d6 d<b ^dt 

ds ds ds ds 



dd 
ds 



0, 



where A, B, C, and V are known functions. By taking initially 6 — 6q — | , 

' d6 s 

0, we get from the above equation 



and 



ds 



dM_ 
ds 2 



= 0, 



then 



ds 



(4.1) 



It is clear from this solution that the motion of a test particle is a planer mo- 
tion. Now restricting ourselves by taking the plane 6 = n/2 for an eclipsing 
binary. In this case the line element will become 



ds 2 = -e^dr 2 - r V# 2 - 2ur 2 e fi d(j>dt + (e v - w 2 rV) dt 2 
For x 3 = 0, the equation of motion becomes, 



d 
ds 



-2r V0 - 2ur z eH 



(4.2) 



(4.3) 



i 9/1 dv 

where 11& = — — , and V& = —. This leads to 



rV (j> + cut) = - l - J e v (^ - /^) idt + Jf, 



(4.4) 



where J2? is an arbitrary constant. 
Similarly, for x 4 = t, we get 



d 
ds 



e v (y t - fit) i 2 



-2wrV0 + 2 (e" - u; W) i 
The above differential equation leads to 

; M (J) + cut) - e v i = ~ J e u [y t - fi t ) idt-S, 
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cur 2 e r 



(4.5) 



(4.6) 



where $ is an arbitrary constant. 



Recalling equations (2.12.1) and (2.12.2), one can show that: 



(4.7) 



Multiplying (|4.4|) by the constant u, combining it with (4.6) and solving the 

(4.8) 



equations for t, we get 

i = {S + wJSf) e~ v . 
Similarly, solving for <j), we get 

Y 2r 

where 



(4.9) 



It is easy to show that 
being, therefore 



2r 2 
O 



(z/0 - nt) dt. 

which can be ignored for the time 

(4.10) 



to {£ + uSf) e~ v + O 



Recalling equation (|4.2), and substituting from (4.8) and (4.10), we get for r 

{£ + uo^f e-**-" - e~» (^B + ^e~^j, (4.11) 
where the parameter B is defined as 



B 



0, for a photon; 

1, for a material particle. 



In what follows, we are going to extract some physical features from the 
above results by considering some special cases. 



4.1 Boundary conditions and physical meaning 
• The field far from the source: 

Assuming that the source of field is too distant from the observer, then we can 
consider that r ~ n ~ r^. Also the total mass of the system is m = mi + m2. 
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By taking the approximation r » m so that O ( — ) — > 0, therefore 



2 

m 



e" = e - " « 1 



2m N 



and 



0=0, (4.12.1) 
r 

* - ^- W 7^y (412 - 3) 



-•■^) 2 7J^ ';,J- (4-12.4) 

It is clear that the above equations are similar to the Schwarzschild case 
except for some additional terms depending on the angular velocity u. Also, 
from our experience in classical mechanics it is clear that the constants of 
integrations 5£ and $ are, respectively, representing the angular momentum 
and the energy of the moving test particle. 

• Static Field: 

Taking u = 0, we get the following set of differential equations: 

= 0, (4.13.1) 

1 « -7 — - — r-, (4.13.2) 
r J 

* « ^» ( 4 - 13 - 3 ) 

<"-(i-^)(b + £). (4.13.4) 



2m \ / _Sf 
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which is now identical to the motion in Schwarzschild field. 

After all, our attempt is to show how some measurable quantities (e.g. 
redshift and time delay in the field of binary systems) can be extracted by 
using of the solution of the equations of motion. 

5 Prospective 

The modified model enables us to go over to radiative spacetimes, considering 
some new developments, especially, impulsive gravitational waves generated 
by boosting various particles. Also, the more useful benefit of this model is 
its capability to calculate the redshift of the pulsations of highly relativistic 
systems (binary pulsars) like the recently discovered double-pulsar system 
PSR J0737-3039A & B [3j[7]. Furthermore, our model enables us to calculate 
the energy and the angular momentum of the binary pulsars by using one of 
the famous methods (e.g. Moller energy-momentum complex). In this way 
we can estimate the decay of the orbital period, coalescence rate and the 
gravitational radiation of the compact binary systems. 
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